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Abstract 

We report a remarkable type of bifurcation: by varying real parameters, unstable 
complex orbits may become stable over wide parameter ranges. Thus, phase dia- 
grams obtained by analizing solely the stability of real solutions may be incomplete. 
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The purpose of this paper is to report a remarkable new type of bifurcation: 
unstable complex orbits may be stabilized by varying real model parameters. 
In other words, by varying real parameters it is possible to stabilize "complex 
phases" in phase- diagrams. This surprising fact is shown for the paradigmatic 
example of a multidimensional dissipative dynamical system, the Henon map 
(x, y) i — ^ (a — x 2 + by,x). The parameter space of the map contains a wide 
domain of real parameters a and b where it is possible to find complex 'ghost' 
solutions which are stable. 

This new bifurcation is of importance in the construction of phase diagrams, 
usually constructed by sweeping real parameters and studying the set of real 
solutions, since domains of complex stable motions might be missing in them. 
Another interesting implication is that the plethora of ghost (complex) orbits, 
so fundamental nowadays in quite different fields [1,2, 3], may be subdivided 
dichotomically into unstable and stable ghosts, pointing to the necessity of in- 
vestigating the effect of complex stability in all physical applications. In atomic 
physics, for instance, the stabilization of complex ghost orbits is expected to 
allow sum re- arrangements in trace formulas [4]. From the exact analytical re- 
sults reported here one can show that the Henon map displays Naimark-Sacker 
bifurcations and, consequently, supports quasiperiodic behaviors[5]. 
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The possibility of stabilizing complex orbits seems not to have been considered 
before [6,7], perhaps because the algebraic varieties involved are of very high 
degrees, exceeding by far those studied by mathematicians [8, 9]. The stabiliza- 
tion of complex orbits was not considered in the classic work of ArnoldflO]. 

As shown recently[ll,12], one may always reduce the equations of motion 
of any algebraic dynamical systems to a pair of polynomial equations: (i) 
P(x, a) = 0, parameterizing simultaneously all orbits of any given period in 
terms of the sum a of orbital points, and (ii) §(<r) = 0, defining the values 
of a as a function of model parameters. The degree of §((x) tells the quantity 
of independent solutions available. In addition, it is also possible to write the 
secular equation ruling the stability of the system as a function of o and of 
model parameters. Following Ref. [11], the polynomials providing complete 
information about all possible period-6 orbits are 



F(x, a)=x 6 - ax 5 + 9 4 (a)x 4 - 9 3 (a)x 3 + 9 2 (a)x 2 - 9 1 (a)x + 9 (a), (1) 
%) = E ©iMK (2) 

where a = X\ + x 2 + £3 + rr 4 + x 5 + x 6 is the sum of the orbital points. 
The coefficients 9i(o~) = 9i(a, b; a) are the standard symmetric functions of the 
roots xe (orbital points). The coefficients 9i(a) and Oj(a, b) are given explicitly 
in the Appendix at the end of the manuscript. Additionally, orbital stability 
is ruled by the following quadratic for the eigenvalues A: 

L(A) = A 2 - [2b 3 + ( M + 2Ui ) b 2 ~2^b + 649 ] \ + b 6 (3) 

where A/5, not contained in Eqs. (l)-(2), is given in the Appendix. Equations 
(l)-(3) result from quite long algebraic manipulations which were performed 
automatically on a computer using specially devised ad hoc routines. 

The degree of Eq. (2) tells that for any set (a, b) of parameters there are nine 
possible period-6 orbits, not all necessarily different. The actual orbits are 
found by substituting the nine roots of Eq. (2) into Eq. (1). Since for real 
parameters all 6j(a, b) are real, (i) there is always at least one real value of a, 
and (ii) complex values of a must always appear in conjugate pairs. 

To get a feeling about the nature of the foliated surface defining a values, 
Fig. (1) shows the real roots of S(er) = as a function of a, for b = —0.98. 
The four points Tj indicate the location of tangent bifurcations. Of special 
interest is the locus or defined by that root of Eq. (2) that is always real. 
Along this locus two different phenomena occur. First, one finds the familiar 
3^6 period-doubling bifurcation, indicated by PD. The doubled orbit is 
stable in a small interval between two vertical dashed lines that is too small to 
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Fig. 1. Real roots of S(cr) = as a function of a, for b = —0.98. Stable complex 
orbits exist between N\ and N 2 on the locus ctr, the locus of the real root that is 
always present. PD refers to the period-doubling, Tj to tangent bifurcations. The 
figures on the right show details hard to see on the left. See text. 

be discernible in the figure. Second, along <jr we find a remarkable new type of 
bifurcation arising from the multivalued character of o~r between the points N\ 
and N 2 . In this interval there are three real roots a, which define three complex 
orbits. The complex orbit defined by the "middle branch" interconnecting Ni 
and N 2 , is stable inside a region resembling a "bow-tie" [see Fig. (3) below] 
precisely where o~r displays a fold [in the interval between Ni and N 2 in 
Fig. (1)] This shows that folds are not always necessarily connected only with 
tangent bifurcations. 

Full lines in Fig. (2) show how the singularities in Fig. (1) evolve when b 
changes. Dotted lines, obtained by investigating eigenvalues, delimit stabil- 
ity domains in the usual way. For reference, Fig. (2) also displays the in- 
terval where period-1 orbits (fixed points) are stable. The new bifurcation 
being reported here occurs inside the box, shown magnified in Fig. (3). As is 
known, tangent bifurcations may be located analytically from the discriminant 
of Eq. (2) with respect to a. 

Figure (3) shows the stability domain of complex orbits along with a line "1" 
marking the border where stable orbits of period-1 are born when a increases. 
Recall that physically meaningful solutions exist only in the interval —1 < 
b < 1. The scenario in the region above the line b = — 1 is divided into four 
domains labeled Ri and displays the following characteristics. All nine period- 
6 orbits are complex in Ri, despite the fact that one of them corresponds 
to a real a. Moving from Ri into R 2 , one finds the complex orbit stabilizing 
bifurcation, when two additional values of o become real (see curve or in 
Fig. (1)) but their corresponding orbits remain complex, one of them being 
stable. In i? 3 there are three orbits associated with real values of a, two orbits 
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Fig. 2. Stability domains and singular- Fig. 3. Blow-up of the box in Fig. 2. 
ities for period-6 motions. See text. Stable complex orbits exists in R 2 . 

being real, one stable and one unstable. The complex orbit is unstable, being 
the same orbit that will give rise to a stable orbit following the 3 — > 6 period- 
doubling, when a increases. When moving from R3 to R4, the stable real 
orbit looses its stability. The intersection of the "1" line with R 2 , located at 
(a*, 6*) ~ (—0.94832, —0.94764), is defined by algebraic numbers of degree 38, 
a quite high degree. This intersection has a duaZ[ll] at 6 = 1/6* ~ —1.0552 and 
a ~ -1.0560. The tip of R 2 , located roughly at (a, 6) = (-0.95449, -0.93348), 
is defined by algebraic numbers of degree 76. 

All in all, the exact expressions reported here, result of long and elaborate 
algebraic computations, reveal the existence of a new sort of bifurcation that 
occurs among complex trajectories, transforming unstable ghost orbits into 
stable complex orbits. Our Eqs (1-3) contain many additional features that 
will be considered in a future publication [5]. 

AE is a doctoral fellow of the Conselho Nacional de Desenvolvimento Cientffico 
e Tecnologico (CNPq), Brazil. JACG is a Research Fellow of the CNPq. This 
work was sponsored by CAPES (Brazil) and GRICES (Portugal). 

Appendix: The Coefficients 9i(a,b;a) and Oj(a, 6) 

The coefficients 9i(a, 6; a) needed to define all period-Q orbits are the following: 
4 = [<7 2 + (l-6)0--6a]/2, 

# 3 = [2a 3 + 6(1 - b)a 2 - 32aa - M(l - b)/V)/12, 

6 2 = (l/24){(x 4 + 6(1 - b)a 3 - [28a - 3(1 - 6) 2 }a 2 - 6(1 - 6) x 

[4a - (1 + b 2 )]a + 36[2a - (1 + 6 2 )] - 2(1 - h)oM x jV + 36A/" 2 /X?}, 
6 1 = (l/120){a 5 + 10(1 - b)a 4 - 5[8a - 3(1 - 6)> 3 - 10(1 - 6)x 

[10a - 3(1 + b 2 )]a 2 + 12[22a 2 - 13a(l + b 2 ) - 26(1 - 6) 2 ]a 
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-(1 - 6) [5a 2 + 5(1 - b)a - \%a\Mi/V 
+3(56(7 - 2(1 - 6 3 )A/" 2 /£> - 6bM 3 /V}, 
6 = (l/720){a 6 + 15(1 - b)a 5 - 5[10a - 9(1 - 6) 2 ]a 4 
-5(1 - b) [16a - 76 2 + 26 - 7]a 3 

+2[272a 2 - 9(276 2 - 106 + 27)a + 9(5 6 2 - 86 + 5)(1 - 6) 2 ]a 2 
+24(1 - 6)(15a - 256 2 + 266 - 25)aa - 360a 2 (2a - 3(1 + 6 2 )) 
-2[(1 - 6)(5a 3 + 15(1 - b)a 2 - Uaa) + 15(6 4 - 6 3 
+66 2 - 6 + l)]Hi/V - 9156a 2 - (1 - 6)(46 2 - 6 + 4)<r - 106a]W 2 /£> 
+3QbaM 3 /V - 3606(1 + b 2 )MjV - 5(1 - 6) 2 (M/^) 2 /2}, 

where the following abbreviations are used: 



V = 3a 4 - 7(1 - 6)a 3 - (12a - 116 2 - 56 - ll)a 2 + (1 - 6) x 
(16a - 196 2 - 346 - 19) <r - (1 - 6) 2 (4a - 86 2 - 136 - 8), 

Mi = a 6 - 2(1 - 6)a 5 + 2(2a + 56 2 - 226 + 5)<7 4 

-4(1 - 6) (18a + 56 2 - 206 + 5)cr 3 - [32a 2 - 8(226 2 - 56 + 22)a 
+3(96 2 + 46 + 9)(6 2 + 46 + l)]a 2 + 2(1 - 6) [64a 2 - 4(296 2 + 466 
+29)a + 516 4 + 1326 3 + 2586 2 + 1326 + 51]<r - 12(1 - 6) 2 x 
[8a 2 - (136 2 + 286 + 13)a + (6 2 + 6 + 1)(86 2 + 136 + 8)], 

M 2 = 3a 6 - 4(1 - b)a 5 - 12(3a - Ab)a A + 2(1 - 6) (48a + 56 2 

-266 + 5)(T 3 + [96a 2 - 24(76 2 + 36 + 7) a + 136 4 + 686 3 + 1266 2 
+686 + 13](7 2 - 2(1 - 6) [64a 2 - 88(1 + 6) 2 a + 276 4 + 486 3 + 906 2 
+486 + 27]^ + 4(1 - 6) 2 (8a - 96 2 - 246 - 9)a, 

M = a 7 - 2(10a - 276)a 5 + 2(1 - 6) (20a - 76 2 - 5 56 - 7)a 4 + [64a 2 
-16(86 2 + 116 + 8)a + (76 2 + 46 + 7)(96 2 + 106 + 9)}a 3 
+2(1 - 6) [16a 2 + 4(296 2 + 536 + 29)a - (1056 4 + 2116 3 + 3 766 2 
+2116 + 105)](T 2 - 4(1 - 6) 2 [48a 2 - (256 2 + 426 + 25)a - (86 2 
+136 + 8)(76 2 + 46 + 7)}a + 24(1 - 6) 3 (4a - 86 2 - 136 - 8)a, 

7V 4 = a 6 - (1 -b)a 5 - (14a + 6 2 - 176 + l)a 4 + (1 - 6) (34a - 36 2 
-436 - 3)a 3 + [40a 2 - 2(296 2 - 6 + 29)a - 46 4 + 276 3 + 86 2 
+276 - A]a 2 - (1 - 6) [96a 2 - 2(696 2 + 1106 + 69)a 
+166 4 + 876 3 + 1786 2 + 876 + 16] <r + 2(1 - 6) 2 x 
[28a 2 - 9(66 2 + 116 + 6)a + 3(6 2 + 6 + 1)(86 2 + 136 + 8)], 

J\f 5 = a s + 2 (-9 a + 2 6 2 + 17 6 + 2)a 6 + 2 (6 - l)(-6 a + 9 6 2 + 29 6 
+9)(7 5 + (72 a 2 - 12 (3 6 2 + 22 6 + 3)a + 186 6 2 + 82 6 + 59 6 4 
+82 6 3 + 59)a 4 + 2 (6 - 1) [8 a 2 - 12 (3 6 2 + 16 6 + 3)a 
+ (6 2 + 6 + 1)(111 6 2 + 116 6 + lll)]cr 3 + 2[-32 a 3 - 24 (2 6 2 
-9 6 + 2)a 2 + (115 6 4 - 56 6 3 - 22 6 6 2 - 56 6 + 115)a + 104 6 6 
+36 6 5 + 22 6 4 + 226 2 + 36 6 + 104]a 2 + 4 (6 - 1) [-64 a 3 
+4 (21 6 2 + 62 6 + 21)a 2 - (11 6 4 + 196 6 3 + 40 2 6 2 + 196 6 + ll)a 
+(6 2 + 6+ l)(16 6 4 + 87 6 3 + 178 6 2 + 8 7 6 + 16)]<r 
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+8(6 - l) 2 [-24a 3 + 5 (11 6 2 + 20 6 + ll)a 2 - 9 (6b 2 + lib 
+6)(6 2 + b + l)a + 3 (8 b 2 + 13 b + 8)(6 2 + b + l) 2 ]. 

The coefficients 0j = 0j(a, b) needed in Eq. 2 to obtain the 9 solutions up = 
a e (a, b) are 9 = 1, 8 = 1 - 6, 7 = -24a + 2(6 2 + 166 + 1), and 

6 = 2(1 - 6) [4a - (76 2 + 126 + 7)], 

5 = 144a 2 - 16(6 2 + 256 + I) a + (496 4 + 526 3 + 2666 2 + 526 + 49), 
4 = -(1 - 6) [112a 2 - 16(6 2 + 276+ l)a 

+ (1756 4 + 3886 3 + 5 1 86 2 + 3886 + 175)], 
3 = -4[64a 3 - 8(6 + 5) (56 + l)a 2 - 2(176 4 - 486 3 - 1726 2 

-486 + 17)a - (7b 2 + 26 + 7)(56 4 + 96 3 - 6 2 + 96 + 5)], 
2 = 4(1 - 6) [64a 3 - 8(156 2 + 386 + 15)a 2 + 6(236 4 + 486 3 + 966 2 

+486 + 23)a - (76 4 + 336 3 + 556 2 + 336 + 7)(76 2 + 26 + 7)], 
B 1 = 8(1 - 6) 2 [32a 3 - 2(196 2 + 346 + 19)a 2 - 2(266 4 + 236 3 + 96 2 

+236 + 26)a + (6 2 + 6 + 1)(76 2 + 26 + 7)(86 2 + 136 + 8)], 
©o = -16a(l - 6) 3 [16a 2 - (376 2 + 626 + 37)a + 3(6 2 + 6 + 1)(86 2 + 136 + 8)]. 
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